Phys. Chem. Earth (B), Vol. 25, No. 5-6, pp. 471-474, 2000 

Pergamon © 2000 Elsevier Science Ltd 
All rights reserved 

1464-1909/00/$ - see front matter 


PII: $1464-1909(00)00047-2 


The Determination of the Effective Radius of Drops in Water Clouds from 
Polarization Measurements 

A. A. Kokhanovsky 

Department of Chemical Engineering and Chemical Technology, Imperial College of Science, Technology and 
Medicine, London SW7 2BY, U.K. 

Institute of Physics, 70 Skarina Avenue, Minsk 220072, Belarus 


Received 25 August 1999; accepted 28 February 2000 


Abstract. A study of the influence of the solar angle, scattering can be used to determine the thermodynamic 
thickness, and microstructure of a cloudy medium on the phase of clouds. Such a determination is based on the fact 
intensity and the degree of polarization of the reflected that the presence or absence of rainbows depends on the 
solar light was performed in the framework of vector shape of particles in a scattering medium [18, 25-27]. A 
radiative transfer theory. Calculations were carried out rainbow does not exist at all for ice clouds which are 
for water clouds for wavelengths of 443, 670 and 865 collections of highly nonspherical particles of different 
nm. The solar angles were varied in the range 0-90 shapes [3, 18, 26]. This feature has already been used in 
degrees. It was found that the rainbow polarization is satellite retrieval algorithms [25]. It could be applied for 
larger for larger particles. The position of the rainbow the discrimination of cloud and snow fields on satellite 
shifts to larger scattering angles with decreasing particle images as well. 

size. The degree of polarization decreases with the optical The main task of this paper is to show that the rainbow 
thickness and reaches an asymptotic value for semi- scattering can be used not only for the detection of the 
infinite clouds at an optical thickness of about 100. A . shape of particles, but also for the determination of the 
new method for the retrieval of optical thickness and the effective radius of particles a „and the optical thickness 


effective radius of water drops at a single wavelength is 
proposed. This method is based on the measurements of 


the intensity and the degree of polarization of the With this in mind extensive calculations of the degree of 


reflected light at the rainbow angle. Such measurements polarization and the reflection function of cloudy media 
are important also for the discrimination of ice clouds and were performed at wavelengths 443, 670 and 865 nm. 


snow fields on satellite images due to the fact that ice 
clouds and snow do not show enhanced scattering at the 


Tt of turbid media. Some results in this direction have 
already been reported by Breon and Goloub [28]. 


These values correspond to channels (with the possibility 
to measure the Stokes vector components) of the 


o 2000 tev, POLDER radiometer on board the Japanese ADEOS-II 
© 2000 Elsevier Science Ltd. All rights reserved. satellite, 

Introduction II. Radiative and polarization characteristics of water 
The problem of the determination of the size of particles clouds 

in turbid media is of importance for different branches of Polarization characteristics of light, reflected from turbid 
modern science and technology, including clouds layers can be obtained in the framework of vector 
microstucture studies from space. It is a well known fact radiative transfer theory [2, 13, 29]. To do so one should 
that water clouds are important regulators of shortwave solve the vector radiative transfer equation with 
fluxes in the Earth atmosphere. Their influence on the corresponding boundary conditions [13, 19-21]. The most 
intensity of solar radiation, coming to the Part was simple and transparent results can be obtained for normal 
investigated in numerous papers and books [1-17]. illumination of 2 layer or for nadir viewing geometries. 
Another important feature of clouds, namely their ability These permit reduction of the number of Stokes 
to change the polarization characteristice of cident parameters to be calculated and analyzed. For instance, 
radiation, was studied in detail as well [18-21]. This the Stokes vector of the reflected solar radiation does not 
knowledge was applied to the solution of the broad depend on azimuth at nadir viewing angles. Moreover, 
spectrum of inverse problems [2, 3, 5, 7-9, 11, 15, 17, 22- the third and forth components of the Stokes vector are 
25]. For instance, it was found [18, 25] that rainbow equal to zero due to the symmetry of the problem. Thus, 
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the reflected radiation is partially linearly polarized with 
degree of polarization P. The electric vector of the 
reflected light is in a plane which is defined by the 
normal to the light-scattering layer and the direction of 
the reflected light (P<0) or in the direction perpendicular 
to this plane (P>0). We found that for most incidence 
angles the value of P is positive and the reflected 
radiation is partially linearly polarized in the direction 
perpendicular to the meridional plane. The degree of 
circular polarization is equal to zero. This approximately 
holds for near nadir measurements as well. 

Our calculations were performed with a vector radiative 
transfer code, based on the doubling method of radiative 
transfer theory [13, 19, 30] at different light incidence 


angles and the near nadir viewing angle &=2°. The 
relative azimuth ø was equal to zero. Some results of 


the calculations are presented in Figs. 1, 2. 
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Fig. 1 The dependence of the degree of polarization of light reflected 


from water clouds on the incidence angle at 4=865nm, 
3=2°,9=0°,a, =6 ym, and T = 5(1),10(2), 30(3), 100(4) . 
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Fig. 2 The same as in Fig.l, but atA = 443 nm, r=10, and a= 4 


(1), 6(2), and 12(3) um . 


Generally speaking, the reflection function and degree of 
polarization of the light reflected from a turbid layer for a 
fixed observation geometry depends on its optical 
thickness , size of particles , and their refractive index. 
The effective radius of the particlesa,,, defined as the 


ratio of the third to the second moment of the particle size 
distribution, is in the range 4 — 16 um and the optical 
thickness r is in the range 5 - 100 for most of water 
clouds. For instance, during the experiment[11] it was 
found that the optical thickness of clouds was in this 
range with the most probable values approximately equal 
to 10. Thus, we performed calculations at the optical 
thickness r= 5, 10, 30, 100 for the values of ay equal 


to 4, 6, 8, 12 and 16 um. The incidence angle was varied 
from zero to 90 degrees. We used the following particle 
size distribution to find the local optical characteristics of 
cloudy media with the Mie theory [19]: 


f(a) = Ba‘ exp(-9a/a,), 


where B is the normalization constant( i f(a)da =1). 


The dependence of the degree of polarization on the 
optical thickness of a layer and the incidence angle 9, is 


presented in Fig. 1. One can see that the polarization 
decreases with the optical thickness. It has a peak, which 
corresponds to the rainbow scattering region. The 
dependence of P(T ) can be used for the determination of 
the optical thickness of clouds from polarization 
measurements. The most sensitive angles are in the range 
of rainbow scattering. It should be pointed out that this 
technique would be of importance only for small and 
moderate values of the optical thickness. The degree of 
polarization only weakly depends on the optical thickness 
for thick layers. It rapidly approaches its asymptotic value 
for a semi-infinite medium. . 
The dependence of the degree of polarization on the size 
of particles is presented in Fig. 2. One can see that larger 
particles are responsible for larger values of the degree of 
polarization inside the rainbow peak. The maximum of 
polarization moves to smaller scattering angles with 
increasing size of the droplets. 


III. Rainbow polarization as a tool for particle sizing 


The dependence of the degree of polarization at the 
rainbow angle, on the value of a, can be used for the 
solution of the inverse problem of cloud optics. In the 
framework of this new retrieval method one should 
measure the values of the reflection function and the 
degree of polarization at the rainbow scattering angle, 
which is defined by the following equation [31]: 
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2 — p? 
0, =m +2arccosl i i —4arcsin aon . 
3 3n? 


This angle depends on the refractive index n of the 
particles, which is different for different wavelengths. 
One can find that the rainbow angle is equal to 138.5, 
137.6, and 137.2 degrees for wavelengths equal to 443, 
670, and 865 nanometers respectively. Thus, the rainbow 


reflection function 


0 0.05 0.1 


shifts away from the backward direction with 
wavelength. 

It was found that the reflection function only weakly 
depends on the size of droplets in visible. This feature is 
attractive for the optical thickness retrieval, as was 
discussed by King [11]. The measured value of the 
degree of polarization for the optical thickness r, found 
from reflection function measurements, can be used to 
estimate the effective radius of the droplets (see Figs. 2, 
3). 
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Fig. 3 The reflection function and degree ofpolarization of light reflected from water clouds at different values of the optical thickness and effective radius 


of droplets for a wavelength of 670nm (nadir observation, 3, = 36°). 


Note that values of 6, obtained from the equation 


_ presented above, are approximate and do not depend on 
the size of the particles. In reality the position of the 
rainbow is shifted and this shift is proportional to the 


2/3 
value of 5 [31,32]. This follows from Fig. 2 as 
a 
Ee 

well. 

To find values of the effective radius and the optical 
thickness one should construct special look-up tables for 
degree of polarization and reflection function of clouds at 


different a, andr and viewing geometries [17]. After 


that, the size of droplets and optical thickness can be 
calculated from the minimum of the following function: 


F=(nR, -InR,)? +(InP, -InP)’, 


where R,,R,,P,,andP are measured (m) and 
calculated(c) reflection functions and degrees of 
polarization respectively. The values of a, andr can be 


found separately at r > 30, because the curves in Fig.3 
are nearly orthogonal for thick clouds. Moreover, the 
reflection function is almost independent of the optical 


thickness as r — ©. Thus, the value of a ef can be found 


from measurements of P(6,) alone in this specific case. 
It is interesting that the maximum of the rainbow shifts to 
smaller solar angles with decreasing particles sizes (see 
Fig. 2), making the degree of polarization at the rainbow 
angle more sensitive to the radii of droplets. This 
improves the retrieval accuracy. 

The method proposed is readily generalized for any 
viewing angle. The viewing angle &, the relative 
azimuth @, and the solar angle 3, are related by the 
following equation [29]: 


cos @, = sin Ysin 3, cos p—cos Fcos 3, 
at the rainbow geometry. At 9 = 0 , one can obtain: 


9, =x-80,. 


r 
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IV. Conclusions 


A method for the determination of the size of water 
droplets in clouds is proposed. The chief shortcoming of 
the retrieval algorithm proposed is that it needs a special 
geometry for the measurements (rainbow scattering). One 


way around this problem is multidirectional measure- . 


ments, which are performed by the POLDER instrument 
[25]. 

The effective radius, obtained with this approach, is not 
representative of the total geometrical depth of the cloud. 
The value retrieved is mostly related to the microphysical 
parameters of the cloud top. This is also the case for the 
particle sizing approach of Breon and Goloub [28], based 
on measurements of light reflection and polarization in 
the region of secondary rainbows. The information 
obtained with the rainbow scattering technique could be 
of value in addition to results obtained with the use of 
spectral signatures of water clouds [8]. 

Note, that the accuracy of the retrieval can be improved 
if one performs spectral measurements of the degree of 
polarization at the rainbow angle. The angular distance 
between secondary and primary rainbows also provides 
additional information on the cloud microstructure. 
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